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Abstract — In cognitive radio, spectrum sensing is a key 
component to detect spectrum holes (i.e., channels not used 
by any primary users). Collaborative spectrum sensing among 
the cognitive radio nodes is expected to improve the ability of 
checking complete spectrum usage states. Unfortunately, due to 
power limitation and channel fading, available channel sensing 
information is far from being sufficient to tell the unoccupied 
channels directly. Aiming at breaking this bottleneck, we apply 
recent matrix completion techniques to greatly reduce the 
sensing information needed. We formulate the collaborative 
sensing problem as a matrix completion subproblem and a 
joint-sparsity reconstruction subproblem. Results of numerical 
simulations that validated the effectiveness and robustness of 
the proposed approach are presented. In particular, in noiseless 
cases, when number of primary user is small, exact detection 
was obtained with no more than 8% of the complete sensing 
information, whilst as number of primary user increases, to 
achieve a detection rate of 95.55%, the required information 
percentage was merely 16.8%. 

I. Introduction 

Cognitive radio [T), J2] has been known as a novel 
paradigm for improving the utilization of the precious natural 
resource - radio spectrum. Spectrum sensing is a key compo- 
nent in cognitive radio for detecting spectrum holes, which 
are the spectrum channels not used by any primary user. 
Since each cognitive radio (CR) node has only limited local 
observation to the whole spectrum due to various constraints, 
collaborations among CR nodes are important for acquiring 
the complete spectrum information. 

Related work J3)-||5) introduces architecture and network 
protocols using either distributed or centralized approaches 
for collaborative sensing in CR networks. In the distributed 
approach, CR nodes exchange information using local com- 
mon channels through distributed coordination. In the cen- 
tralized approach, a central control entity, e.g. a fusion 
center, gathers sensing information from all the CR nodes 
within a network through the common control channel. In 
both cases, subject to energy constraints, each CR user has 
very limited access to the whole spectrum. Furthermore, 
due to path loss, channel fading, and/or shadowing effects, 
the highly incomplete sensing information transmitted by 
CR users is affected by transmission loss or even error. 
Therefore, the challenge is for the system to obtain complete 
channel states from their incomplete measurements. 



Contributions: 

We assume that a few of the n channels are occupied by 
primary users and channel fading is not known. A system 
model is introduced in which each CR, instead of scanning 
each channel, takes a small number of measurements that 
linearly combine multiple channels, which are sent to the 
fusion center. We allow some transmissions to fail so the 
fusion center receives a subset of the measurements. To 
assess the channel occupancy, we first argue that the matrix 
of all the measurements has a low-rank so it can be recovered 
by solving a matrix completion problem. Next, given the 
recovered matrix, the occupied channels and their fading val- 
ues are reconstructed by solving a joint-sparsity reconstruc- 
tion problem. In our simulations we used FPCA, a matrix 
completion algorithm from iTTOl . and a novel joint-sparsity 
algorithm. In particular, in noiseless cases, when number of 
primary user is small, exact detection was obtained with no 
more than 8% of the complete sensing information, whilst as 
number of primary user increases, to achieve a detection rate 
of 95.55%, the required information percentage was merely 
16.8%. 

This paper is organized as follows: In Section HH the 
system model is given. Algorithms for collaborative sensing, 
including the matrix completion and joint-sparsity parts, are 
described in Section [HI] Simulation results are presented in 
Section ITVl and conclusions are drawn in Section IVl 

II. System Model 

Suppose there are m CR nodes locally monitoring a subset 
of n channels where m < n. A channel is either occupied by 
a primary user or unoccupied, corresponding to the states 1 
and 0, respectively. Assume that the number s of occupied 
channels is much smaller than n. Our task is to tell the 
occupied channels from the CR nodes' observations. 

The proposed approach for collaborative sensing is de- 
picted in Figure Q] Instead of sensing one channel at a time, 
each CR takes measurements of multiple channels using the 
equipped frequency selective filters, and the measurements 
are sent to the fusion center. Supposing that totally p reports 
regarding the n channels are sent from a CR to the fusion 
center, we model this process by a p x n filter coefficient 
matrix F. Let an?ixn diagonal matrix R present the states 
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Fig. 1. Collaborative Sensing using Matrix Completion 



of all the channels using and 1 diagonal entries, indicating 
the unoccupied and occupied states, respectively. There are 
s entries of 1 in diag(R). In addition, channel gains that 
affecting the CR nodes are described in an m x n channel 
gain matrix G given by 



G, 



(dij)~ a/2 \hiA 



(1) 



where d^j is the distance from the primary transmitter using 
j th channel from the i th CR node, a is the propagation 
loss factor, and hi j is the channel fading. The measurement 
reports sent to the fusion center can be written as a p x to 
matrix 
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Measurement matrix M has the following two important 
properties (5) required for completion from partial entries: 

1) Low Rank: rank(M) equals to s, which is the number 
of prime users in the network and is usually very small. 

2) Incoherent Property: Generate F randomly (subject 
to hardware limitation). From (Q~|) and the fact that R 
has only s nonzeros on the diagonal, M's SVD factors 
U, S, and V satisfy the incoherence condition liTTl 

• There exists a constant /in > such that for all 

i G \p\, j € [m], we have Y?k=i U ?,fc - Mos, 
Efc=l Y lk < Vos- 

• There exists fix such that | Yfk=i ^» fc^feVj k |< 
Mis 1/2 - 

M is in general incomplete because of transmission 
failure. Moreover, each CR might only be able to collect 
a random (up to p) number of reports due to the hardware 
limitation. Therefore, the fusion certain receives a subset 
set EC [p] x [m] of M's entries. We assume that the 



received entries are uniformly distributed with high prob- 
ability. Hence, we work with a model in which each entry 
shows up in E identically and independently with probability 
e/y/p x to. Given E pxm , the partial observation of M is 
defined as a p x to matrix given by 
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if (i,j)eE, 

otherwise. 



(3) 



We shall first recover the unobserved elements of M from 
M £ . Then, we reconstruct (RG T ) from the given F and M 
using the fact that all but s rows of (RG T ) are zero. These 
nonzero rows correspond to the occupied channels. Since p 
and to are much smaller than n, our approach requires a 
much less amount of sensing and transmission, compared 
to traditional spectrum sensing in which each channel is 
monitored separatively. 

III. CR Sensing Matrix Completion Algorithm 

In previous research on matrix completion i9l- lfT2l . it has 
been proved that under some suitable conditions, a low-rank 
matrix can be recovered from a random, yet small subset of 
its entries by nuclear norm minimization: 
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where |jM||* denotes the nuclear norm of matrix M and r is 
a parameter discussed in Section UlI-BI below. For notational 
simplicity, we introduce the linear operator V that selects 
the components E out of a p x n matrix and form them into 
a vector such that IIT^M 



VM E \\l 



Mf^ l 2 . The adjoint of V is denoted by V* . 

For our problem, we adopt FPCA by Ma et al. in iflOll , 
which appears to run very well for the relatively small- 
dimensional application we focus on. In the following sub- 
sections, we describe this algorithm and the steps we take for 
nuclear norm minimization. We further discuss the stopping 
criteria for iterations to acquire optimal recovery. Finally we 
show how to obtain RG T from M. 

A. Nuclear Norm Min. via Fixed Point Iterative Algorithm 
FPCA is based on the following fixed-point iteration 

- Y k = M k - S k V*(VM k - VM E ) 



M k+1 = S TSk (Y k ) : 



(5) 



where 6k is step size and S a (-) is the matrix shrinkage 
operator defined as follows: 

Definition 1: Matrix Shrinkage Operator S a (-): As- 
sume M G W xm and its SVD is given by M = 
Udiag(cr)V T , where U G R pxr , a G W + , and V G R mxr . 
Given a > 0, S a (-) is defined as 

S T (M) : =Udiag( SQ ( ( 7))V T (6) 

with the vector s Q (er) defined as: 

s a (x) := max{x — a, 0}, component-wise. (7) 



Simply speaking, SV(M) reduces every singular values 
(which is nonnegative) of M by r; if one is small than a, 
it is reduced to zero. 

To understand ©, observe that the first step of (0 is 
a gradient-descent applied to the second term in (01 and 
thus reduce its value. Because the previous gradient-descent 
generally increases the nuclear norm, the second step of (0 
reduces the nuclear norm of Y . Iterations based on (0 
converge when the step sizes 5^ are properly chosen (e.g., 
less than 2 or by line search) so that the first step of (f5]) is 
not "expansive" (the other step is always non-expansive). 

The second step of © requires computing the SVD 
decomposition of Y k , which is the main computational cost 
of ([5]). However, if one can predetermine the rank of the 
matrix M, or have the knowledge of the approximate range 
of its rank, a full SVD can simplified to computing only a 
rank-r approximation to Y fc . Specifically, the approximate 
SVD is computed by a fast Monte Carlo algorithm developed 
by Drineas et al. lfl3l . For a given matrix A S I mx ™ and 
parameters k s , this algorithm returns an approximation to the 
largest k s singular values corresponding left singular vectors 
of the matrix A in a linear time. 

B. Stopping Criteria for Iterations 

We tuned the parameters in FPCA for a better overall per- 
formance. Continuation is adopted by FPCA, which solves a 
sequence of instances of easy to difficult, corresponding 
to a sequence of varying (large to small) values of r. The 
final r is the given one but solving the easier instances of 
gives intermediate solutions that warm start the more 
difficult ones so that the entire solution time is reduced. 
Solving each instance of requires proper stopping. We 
use the criterion: 

ii M fc+i _ M fch 



< mtol 



(8) 



max{l,||M fc || F } 

where mtol is a small positive scalar. Experiments shows 
that le~ 5 is good enough for obtaining correct detections. 

C. Channel Availability Estimation Based on the Complete 
Measurement Matrix 

Since F has more columns than rows, directly solving 
X := RG T in (Hi from given M is underdetermined. 
However, each row Xi of X corresponds to the occupancy 
status of channel i. Ignoring noise in M for now, Xi contains 
a positive entry if and only if channel i is used. Hence, most 
rows of X are completely zero, so every column X.j of 
X is sparse and all X. j's are jointly sparse. Such sparsity 
allows us to reconstruct X from ([1) and identify the occupied 
channels, which are the nonzero rows of X. 

Since the channel fading decays fast, the entries of X 
have a large dynamic range, which none of the existing 
algorithms can deal with well enough. Hence, we devel- 
oped a novel joint-sparsity algorithm briefly described as 
follows. The algorithm is much faster than matrix completion 



and typically needs 1-5 iterations. At each iteration, every 
column X.j of X is independently reconstructed using the 
model niin{^i I : FX.j = M.j}, where M.j is 

the jth column of M. For noisy M, we instead use the 
constraint \\FX.j — M.j\\ < a. The same set of weights u>i 
are shared by all j at each iteration, wt is set to 1 uniformly 
at iteration 1. After channel i is detected in an iteration, Wi 
is set to 0. Through u>i, joint sparsity information is passed 
to all j. Channel detection is performed on the reconstructed 
I.j's at each iteration. It is possible that some reconstructed 
X.j is wrong, so we let larger and sparser I.j's have more 
say. If there is a relatively large Xij in a sparse X.j, then 
i is detected. We found this algorithm to be very reliable. 
The detection accuracy is determined by the accuracy of M 
provided. 

IV. Simulation 

According to FCC and Defense Advance Research 
Projects Agency (DARPA) reports OH, (HI data, we chose 
to test the proposed algorithm for spectrum utilization effi- 
ciency with such settings: at certain times, the number of 
active primary users varies from 1 to 3 on a given set of 
100 channels with 20 CR nodes collaboratively detecting 
the occupied channels. The tests were performed at different 
sampling rates with different numbers of prime users (i.e., 
occupied channels). We define sampling rate as 

No. received measurements at the fusion center 
No. channels x No. CRs 
where (No. channel x No. CR) is the amount of total 
sensing workload in traditional spectrum sensing. Perfor- 
mance was given in terms of probability of detection (POD) 
according to the definition in information-theory. 

POD = No. Hit/ '(No. Hit + No. Miss); 

where No. Hit is the number of successful detection of the 
appearance of primary user(s), while No. Miss is the number 
of miss detections of the appearance of primary user(s). 

Figure [2] shows the probability of detection at different 
sampling rates when there is no noise. In the relatively 
low spectrum usages cases, i.e. one out of 100 channels is 
occupied by a primary user, the proposed scheme enabled 
100% detection probability at a sampling rate as low as 8%. 
With 3 primary users appearing in the network, the detection 
task becomes harder. However, the algorithm still managed 
to realize a detection probability of higher than 95% at a 
low sampling rate of 16.8%. 

Figure [3] shows the probability of detection at different 
sampling rates when the received signal is corrupted by 
Gaussian noise with a signal to noise ratio of 35dB. Col- 
laborative detection becomes even harder, we need higher 
sampling rate for exact primary user detection. Simulation 
result shows that, with 14.4% sampling rate, 1 or 2 primary 
user(s) can be detected exactly. As the number of primary 
users increases, more samples are needed. 
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Fig. 2. POD vs. Sampling Rate (Noise-Free) 

POD vs. Sampling Rate (SNR=35dB) 
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Fig. 3. POD vs. Sampling Rate (SNR=35dB) 

As we can see from the simulation results, the proposed 
approach achieves very high probability of detection at 
an extremely low sampling rate compared to traditional 
spectrum sensing. We plan to improve our algorithms and 
perform broader experiments. The results will be reported in 
a forthcoming journal paper. 

V. Conclusions 

In this paper, we propose a collaborative spectrum sensing 
approach to detect spectrum holes in a cognitive radio 
network. We model the collaborative detection problem as a 
matrix completion problem in which partial observations of a 
matrix enable its faithful reconstruction. We solve the matrix 
completion problem by the recent algorithm FPCA and 
estimate channel availability based on joint sparsity recovery. 
Performance of the proposed approach was tested for fre- 
quency utilization efficiency ranging from 1% to 3%. In the 
noiseless cases, exact detection was obtained with no more 
than 8% of the complete sensing information, whilst as the 
number of primary user increases, to achieve a detection rate 



of 95.55%, the required information percentage was merely 
16.8%. In the noisy cases (SNR: 35dB), less than 15% 
samples enabled exact detection of small numbers of primary 
users. To summarize, the proposed approach significantly 
reduces the amount of sensing and transmission workload 
of cognitive radios for wide range spectrum sensing. 
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